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Abstract 

Due to effective enhancement of the Coulomb coupling strength in the vacuum-suspended 
graphene, the system may turn from a semimetal into an insulator by the formation of 
a gap in the fermionic spectrum. This phenomenon is analogous to the spontaneous 
breaking of chiral symmetry in the strong-coupling relativistic field theories. We study 
this "chiral symmetry breaking" and associated collective excitations on graphene in 
the strong coupling regime by taking U(l) lattice gauge theory as an effective model 
for graphene. Both compact and non-compact formulations of the U(l) gauge action 
show chiral symmetry breaking with equal magnitude of the chiral condensate (exciton 
condensate) in the strong coupling limit, while they start to deviate from the next-to- 
leading order in the strong coupling expansion. Phase and amplitude fluctuations of 
the order parameter are also investigated: in particular, a mass formula for the pseudo- 
Nambu-Goldstone mode (7r-exciton) , which is analogous to Gell-Mann-Oakes-Renner 
relation for the pion in quantum chromodynamics (QCD), is derived from the axial Ward- 
Takahashi identity. To check the applicability of the effective field theory description, 
typical energy scales of fermionic and bosonic excitations are estimated by identifying 
the lattice spacing of the U(f ) gauge theory with that of the original honeycomb lattice 
of graphene. 



1. Introduction 

Graphene, a monoatomic layer material of carbon atoms with honeycomb lattice 
structure, has attracted a great interest both in theoretical and experimental sides. One 
of its important features is the "Dirac cone" structure: since the charge carriers on 
graphene at low energy obey the linear dispersion relation around two independent "Dirac 
points" in the 1st Brillouin zone they can be described as Dirac quasiparticles with 
U(4) chiral symmetry, which means that their bare mass m = [2|. Thus graphene 
shows a gapless spectrum and the density of states vanishes at zero-energy, which means 
that monolayer graphene behaves as a 2-dimensional semimetallic material. It exhibits 
many interesting phenomena: high carrier mobility, anomalous quantum Hall effect, and 
so on Q. 

There is one critical difference between the quasiparticles on graphene and ordinary 
relativistic fermions: the Fermi velocity v F of the quasiparticles on graphene, which is de- 
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termined by physical quantities characteristic to the material, is about 300 times smaller 
than the speed of light c. On the other hand, the propagation speed of the electromag- 
netic field still remains c in the vacuum. Such discrepancy between two velocities leads 
to an effective enhancement of the Coulomb coupling strength, suggesting the possibil- 
ity that the chiral symmetry breaking (exciton condensation) would spontaneously occur 
and that the spectrum would be gapped in the vacuum-suspended graphene [3| ■ Since the 
effective coupling strength is larger than that in QED about 300 times, perturbative ap- 
proach is unsuitable for this system. Thus, we need various non-perturbative treatments, 
which is analogous to the situation in quantum chromodynamics (QCD) [Bj]. Although 
experimental investigation of vacuum-suspended graphene started only recently Q , this 
problem is becoming important in the industrial aspect as well as in the theoretical side. 

There have been various theoretical studies about this problem. Studies with the gap 
equation (Schwinger-Dyson equation) based on 1/r Coulomb potential, first performed 
in 200f, have found that graphene would be gapped in the strong-coupling and small- 
flavor region. They have also estimated critical values of the coupling strength a and the 
flavor number TV [7, 8, 9] . The infrared behavior of the long-range Coulomb interaction is 
investigated by exact renormalization group methods and critical exponents of the phase 
transition are estimated The above studies are focused on the critical behavior of 
the system around the critical values a c and Nc- 

The behavior of the system in the large- N limit has also been studied by 1 /N expan- 
sion. It has been seen that the long-range Coulomb interaction becomes irrelevant and 
that the spectrum would not be gapped at large N, even though the coupling strength 
is large enough [H, 12, 13j |. 

Recently lattice Monte Carlo simulations have been performed based on the Thirring- 
like 4-fermi interaction [l4| and on the U(l) lattice gauge theory 15, 16, 17 1. They have 
calculated chiral condensate (the order parameter for chiral symmetry breaking) with 
varying coupling strength and have estimated the critical coupling strength, with compact 
and non-compact gauge formulations. They have also estimated critical exponents of the 
phase transition and the equation of state in the symmetry-broken region. Due to the 
finite size effect of the lattice, the above lattice Monte Carlo simulations are limited 
to the finite coupling strength (a < oo) and finite bare mass (m > 0) region. The 
behavior of the system in the chiral limit (m = 0), such as the critical coupling value, is 
estimated by extrapolation with the help of the equation of state, introduced in analogy 
with (3+l)-dimensional QED. 

In this paper, we investigate the behavior of this system in/around the strong cou- 
pling (a — > oo) and chiral (m = 0) limit analytically by strong coupling expansion of 
the square-lattice-regularized effective gauge theory model. This method is one of the 
suitable methods to investigate the non-perturbative behavior of gauge theories, such as 
QCD IH lil 20j j , and has been first applied to the graphene effective model in our pre- 
vious paper [21]. U(l) lattice gauge theory is employed, as a low-energy effective model 
for monolayer graphene with Coulomb interaction at zero temperature. An expansion 
around the strong coupling and the chiral limit is performed, which is referred to as 
"strong coupling expansion" of lattice gauge theory. Then we calculate the value of the 
exciton condensate up to the next-to-leading order (NLO) in the strong coupling expan- 
sion, and show the spontaneous chiral symmetry breaking in/around the strong coupling 
limit. Results from the compact and non-compact formulations of the gauge action are 
compared. These results coincide in the strong coupling limit due to the suppression 
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of the pure gauge term, and deviate from the NLO in the strong coupling expansion. 
Our analytic results are also compared with the results from numerical simulations in 



We also investigate two fluctuation modes (bosonic collective excitations) of the exci- 
ton condensate, which have not been studied in the preceding literatures. One of them, 
which we refer to as "7r-exciton," behaves as a pseudo-Nambu-Goldstone boson emerging 
from the spontaneous symmetry breaking. A mass formula for the 7r-exciton is derived 
from axial Ward-Takahashi identity, which is analogous to the Gell-Mann-Oakes-Renncr 
(GMOR) relation for pions in QCD [23 |. The mass of the other one, the "cr-exciton," is 
also calculated, up to the NLO in the strong coupling expansion. In order to check the 
applicability of the Dirac fermion description for graphene at low energy, we compare the 
excitation energies (masses) of the above fermionic and bosonic excitations to the typical 
energy scale from the interatomic spacing of the original honeycomb lattice of graphene. 

This paper is organized as follows: in Section [21 we briefly review the continuum 
effective model and square lattice model for monolayer graphene. In Section|31 we perform 
strong coupling expansion with the lattice effective model and derive the free energy of 
this system up to the NLO, both with the compact and non-compact formulations of 
the gauge action. We calculate the value of the chiral condensate from this effective 
potential. In Section 21 ir- and cr-exciton modes are investigated. Masses of these modes 
are calculated, and a GMOR-like mass formula is derived. Finally in Section \5[ we 
summarize our work and present several problems for future studies. 

2. Low-energy effective model 

In this section, we briefly review the effective field theory which can well describe the 
low-energy electronic behavior of monolayer graphene. We start from the tight-binding 



Here, a\ a and b\ b denote creation and annihilation operators of electrons on the trian- 
gular (Bravais) sublattices A and B respectively, which constitute the honeycomb lattice 
of graphene. r takes the positions on the A sublattice. (i = 1, 2, 3) denote the relative 
position of one B-site from its neighboring A-site, and t corresponds to the hopping am- 
plitude between two neighboring sites. This Hamiltonian has vanishing eigenvalues at 
two independent points K±, which are called "Dirac points," in the momentum space [l|. 
Since the dispersion relation can be linearized around the Dirac points, electron/hole ex- 
citations on graphene at low energy can be described as massless Dirac fermions [2| . Since 
there are 8 degrees of freedom corresponding to 2 (number of sublattices) x 2 (number of 
Dirac points) x 2 (up and down spin) , the quasiparticles on monolayer graphene can be 
formulated by two "flavors" of four-component Dirac fermions. Thus we can construct 
"Dirac spinors" around two Dirac points K± as 



Ref.[l2|. 



Hamiltonian [l[ 




(1) 



reAi=l,2,3 



Vv(p) 



/ a(K++p) \ 
6(K+ + p) 
6(K_+p) 
V a(K_+p) J 



(2) 
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Here we take |p| much smaller than |K±|, which we refer to as "low-energy approxima- 
tion." g =r, I denotes the original spin of the quasiparticle, which is treated as "flavor" 
degree of freedom of the 4-component spinor. 

2.1. Effective model in continuum limit 

In order to incorporate the Coulomb interaction between the quasiparticles, we add 
the electromagnetic field, namely U(l) gauge field, in the effective action. The model is 
described by so-called "braneworld" or "reduced QED" [23l|-like model action 0, 12 1 



] f dx^ {p f (D[A] + m)^ f + ^ 1 Yl fit (1) (U-Hfi (3) 

** & ii i, — 1 o ^ a .** 



E 

/ " J "=1,2, 3,4 ' 



in the Euclidean space-time, where the natural unit (h = c = 1) is taken. In this model, 
the fermion quasiparticles run in the (2+l)-dimensional plane x^ = (t(= X4), x\, X2) 
with the Fermi velocity v F , while the U(l) gauge field propagates in the (3+l)-dimensional 
space = (t(= X4), x\,X2, £3) with the speed of light c(= 1). The Fermi velocity is de- 
fined by two parameters characteristic to graphene lattice, v F — (3/2)£a Hc = 3.02 x 10~ 3 , 
with the interatomic spacin g of the honeycomb lattice a Hc = |sj| = 1.42 A and the hop- 
ping amplitude t ~ 2.8 eV [24|. The above value of v F is the physical value observed 
on Si02 substrate. The "flavor" / of the 4-component Dirac spinor tpf runs from 1 to 
the "number of flavors" N. In this paper we specifically focus on N — 2 case, which 
corresponds to the spin degree of freedom of the monolayer graphene. 
The Dirac operator in Eq.© is defined as 

D[A] = 74(^4 + iA A ) + v F J2 7;(9; +iAi), (4) 

i=l,2 

where denotes /Lt-th component of the U(l) gauge field. Since the fermions are confined 
in the layer, z-componcnt does not appear in the Dirac operator. The gauge coupling 
constant for the vacuum-suspended graphene g 2 — e 2 /eo, with e being the electric charge 
and eo being the dielectric constant of vacuum. If the layer is placed on a substrate, 
interaction strength g 2 is screened by the factor 2/(1 + e), with e being the dielectric 
constant of the substrate [l2j]. Strictly speaking, vp in Eq.Q is a bare value of the 
Fermi velocity which receives a finite renormalization due to interactions. In the strong 
coupling expansion which we employ, bare value of vp is assumed to be as small as the 
physical value. 

The Hermitian 7 matrices obey the well-known Clifford algebra {7^, 7„ } = 25^, and 
the "chirality" matrix 75 is defined by 75 = 74717273- By linearizing the tight-binding 
Hamiltonian in Eq.([T]), three of the gamma matrices 



71= ( V *r 2 J' *=[ in J' *=[ -T 3 ) (5) 

are uniquely obtained corresponding to the Dirac spinor representation in Eq. © . T\ (i = 
1,2,3) are the Pauli matrices which acts on the sublattice (A,B) subspace. Due to the 
absence of z-direction in the tight-binding Hamiltonian, there is a degree of freedom in 
choosing a matrix representation of 73 and 75. One possible choice is 

(6) 





' 


-in \ 




f 




73 = 1 


v *T3 


J 


, 75 = 1 


V T 3 






The bare mass in Eq.Q corresponds to half of the explicit bandgap which originates 
from the difference between on-site energies on A- and B-sublattices. It is suggested that 



such a gap may be formed artificially on epitaxially grown graphene on substrate 25] or 
on graphene nanoribbon and nanomesh [26]. Although it still remains a great question 
in graphene physics how to open a finite spectral gap, in this paper we leave the mass 
term for convenience of calculation. 

2.2. Effective coupling strength and instantaneous approximation 

To set the Dirac operator independent of the Fermi velocity v F , we perform a scale 
transformation in the temporal direction 

t^-t/v f , A 4 -^u f A 4 , (7) 

yielding the effective mass m* = m/v F . This scale transformation renders the Dirac 
operator Lorentz invariant, 

D.[A] = 7 4 (<9 4 + iA A ) + 7i(<9i + iA), (8) 
i=l,2 

and the pure gauge action (the second line in Eq.Q) Lorentz non- invariant, 



Sc 



Y, I dx^idiAo-W (9) 



y i=l,2,3 



+jir E [dx^id^-W. (io) 

Such Lorentz non-invariance leads to the difference between the coupling strength of 
electric (temporal) component and that of the magnetic (spatial) component of the gauge 
field. As for the magnetic component, it is weakened as = v F g 2 , rendering the 
coefficient of Eq. (fTU| sufficiently large to apply saddle point approximation over the 
spatial components of the gauge field. The saddle point gives the vanishing magnetic 
field, B = V x A = 0, and the arbitrariness in choosing A can be absorbed in the 
local U(l) gauge invariance. Therefore we can neglect the spatial components of the 
gauge field Ai (i = 1,2,3). This approximation is usually referred to as "instantaneous 
approximation," because it omits the retardation (magnetic) effect of the electric field. 
Fluctuation correction beyond the saddle point approximation corresponds to the weak 
coupling expansion by g m . In this paper, we work on the saddle point and consider only 
the temporal component of the gauge field. 

On the other hand, the coupling strength of the electric component is enhanced as 
fl* = 9 1 l v F i which is about 300 times larger than the usual Coulomb coupling strength 
in QED. This means that perturbative expansion by as for usual QED does not work 
well in this graphene model, while expansion by the inverse coupling strength 



can be well performed in this regime. From Eg. dill) , we obtain the value of /3 in the 
vacuum 0.0369, while that on the Si02 substrate 0.101. Although the physical value of 
/3 can be shifted due to the renormalization of v F , here we employ the bare value of j3. 
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2.3. Regularization on square lattice 

In order to regularize this theory, we introduce a square lattice so that we can compare 



our results to those of the Monte Carlo simulations Ha, 16, 17 1. We set the lattice spacing 
a comparable to the honeycomb lattice spacing a Hc to reproduce the physical momentum 
cutoff on the honeycomb lattice. Naively the fermionic term of Eq.([3]) is discretized on 
the square lattice with 4-component Dirac spinor tp as 



Sf= ^f(x)[D[U](x,y) 



i>f(y), 



(12) 



where the lattice Dirac operator reads 

D[U](x, y) = - [ S V,*+fi U M( X ) ~ s x, y +p,U*(y)] 7 M - 



(13) 



Here we take the lattice unit, in which all the dimensionful quantities are scaled by a. fx 
denotes the unit vector in rr^-direction. The U(l) gauge field is represented by the link 
variable U fl (x) which corresponds to the link between x and x + fx. The time-like link 
is defined as U±{x) — exp [iO(x)} (— it < 9 < tt), while the space-like links Ui^,3(x) are 
set to unity as a result of instantaneous approximation. Eq. (|12p is the lattice- regularized 
version of the continuum action Eq.Q with the finite lattice spacing a. However, we 
can see in momentum space that such a lattice fermion has 7 unnecessary poles other 
than p M = 0, which means that the lattice action describes 2 3 (= 8) species of fermions 
27 1 . Such species are called "doublers" , and this problem is referred to as the "doubling 
problem." 

In order to avoid such doubling problem, one possible solution is to consider the 
doublers as "real" degrees of freedom. Since the Dirac operator Ea. (fT5]) is diagonalized 
by the transformation 

X (x) = T | 4 7f 7fW, (14) 

the lattice action is written with the single-component Grassmann number xi x )- We 
identify 8 doublers of \ as 4(spinor) x 2 (flavor) deg rees of freedom of the fermion, which 
is called "staggered fermion" formulation [28[ |29(. As a result, the lattice action for 
fermions on graphene reads [151 ] 



\ E (V+(x)-V-{x))+ msf M{x) 



/i=l,2,4 



with fermionic bilinears 



S 
S 

V^(x) = ^r)^(x)x s (x + p,)Ul(x)xs(x), 



(15) 



(16) 



(17) 



where [i = 1,2,4, and the flavor index s of the staggered fermion runs from 1 to N/2. 
In the case of the monolayer graphene, we need no staggered flavor index s, because 
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N/2 = 1 [lj, [l5[ . The staggered phase factors ?7 M corresponding to the Dirac 7-matrices 
are 774(2;) = 1,171(2;) = (-l) r , 772(2;) = (-l) T+Xl . e(x) = (-TT+*i+*2 corresponds to the 
chirality 75. 

As for the pure gauge action term, there are two ways for discretization on the lattice. 
Naively the discretized gauge action is written with the compact link variables U/±{x) as 



c(C) 



E [l-Re(u4x)Ui(x + j))]. (18) 



g * 3=1,2,3 



However, it is known in the normal QED that such compact gauge formulation suf- 
fers from the so-called "monopole condensation" problem, leading to anomalous phase 
transition [301 ] . In order to avoid this problem, we can employ the non-compact gauge 
formulation for this system 15]: 

4 NC) = iE E [«(*)- *(*+i)] a . (19) 

y * xW 3=1,2,3 

We will compare the results from these two formulations later. 

2.4. Chiral symmetry of the system 

In the chiral limit (m — ► 0), the continuum action Eq.Q is invariant under U(4) 
chiral transformation generated by 16 generators {1,<?} ® {1,73,75,737s} for each fla- 
vor, where cr,- (i = 1, 2, 3) are the Pauli matrices acting on the spin subspace. Note that 
invariance under continuous chiral symmetry within the low-energy approximation origi- 
nates from the discrete Z2 symmetry between two triangular (Bravais) sublattices (A,B) 
of the original honeycomb lattice of graphene. If higher order terms in momentum are 
introduced perturbatively, continuous approximate symmetry is eventually broken into 
discrete symmetry. Lack of 73 in the Dirac operator, which comes from the fact that 
the fermion is confined in the (2 + l)-dimensional layer, extends the well-known chiral 
symmetry generated by {1, 75} into such a large symmetry. If the chiral condensate {i^ip) 
obtains a finite expectation value, the chiral symmetry is spontaneously broken as 

U(4) -> U(2) x U(2). (20) 

As for the lattice- regularized model with staggered fermions (Eq. (IT5ll ). the symmetry 
of the system is reduced into U(l) v x U(1) A in the chiral limit: 

U(1) V : {x{x),x{x))^(e^ X (x),e-^x{x)) (21) 
U(1) A : (X(x), X(x)) -> (e^ £ W X W,e^ £ (^W) (22) 

These are remnants of global U(4) chiral symmetry of Eq.Q [31]. Part of the chiral 
symmetry (flavor symmetry) is not conserved in the staggered fermion formulation, which 
is known as "taste breaking." Here we discuss only the spontaneous breaking of the 
U(1) A symmetry. Under the U(1) A rotation, the fermion bilinears transform as M(x) — ► 

(x) , so that the chiral condensate (xx) serves as an order 
parameter for the spontaneous symmetry breaking, 

U(l) v x U(1) A -► U(l) v . (23) 

If (xx) acquires a finite expectation value, the global U(l) v x U(1) A symmetry is broken 
to the mere U(l) v symmetry. 
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3. Strong Coupling Expansion 



Strong coupling expansion on lattice, which employs the inverse of the coupling con- 
stant as an expansion parameter, has been used for analyzing the non-perturbative be- 
havior of gauge theories [l8[ • It was first applied for the pure Yang-Mills theory to show 
the mechanism of confinement of quarks in the strong coupling limit 32|. It was also 



applied to SU(iV) lattice gauge theory with fermions to show the spontaneous breaking 
of chiral symmetry and corresponding meson mass spectra in the strong coupling limit 
successfully (HI- Recently strong coupling limit study of lattice QCD is extended to the 
finite temperature and density region, to map the phase diagram of QCD matter under 
extreme conditions [I^j3|j3i. Finite coupling effect (NLO and next-to-NLO in strong 



coupling expansion) 2(J iff 23 an< ^ the confinement (Polyakov loop) effect [38j are also 
being considered by strong coupling expansion. 

Since QCD is an asymptotic free gauge theory, in which the bare coupling approaches 
zero in the continuum limit, the strong coupling expansion does not give results more 
than a qualitative estimate. On the other hand, in the case of graphene, the lattice 
constant a ~ a Hc is fixed to be a finite value and the bare coupling strength at that scale 
is large enough. Therefore, the strong coupling expansion may provide us with not only 
qualitative estimate but also some quantitative calculations of physical observables. 



3.1. Link integration 

The expansion parameter for strong coupling expansion is defined in Eq. pTj) as f3 = 
l/s*> so that the strong coupling limit corresponds to /3 0. Then, the partition 
function can be expanded by Sq ~ 0(0) around the strong coupling limit as 



Z = J [d X dx}[ 



{-i>G) -S B 



n=0 



= / [dxdxY 



(24) 



Since the integrand of the above equation is written as a polynomial of link variables by 
the expansion, 6 integration can be analytically performed order by order in f3 [l8[ . When 
the link variable e %e and its complex conjugate e~ 10 cancel with each other, the fermion 
self-interaction terms are induced; the terms in which link variables remain vanish away 
through the link integration. 

In the LO (O(/3 )) of strong coupling expansion, there are no difference between the 
compact and non-compact gauge formulations, because the pure gauge term Sq vanishes 
at P = 0. Thus only the fermionic term Sf contributes to the LO effective action S% . 
Due to the Grassmann nature of x, we can rewrite the temporal link integration part as 



(25) 



[do] n 



i-\xs(xy e(x) X s(x + i) 



l + \x s (x + A)e-^Xs{x) 



By integrating out the U(l) compact link variables, we obtain 
exp \-S^ 



M(x+4) CC0) 



M(x) ^ 



,-(x + 4) <f) J|) 



(a) 



(b) 



Figure 1: Induced four-fermion interaction in the strong coupling expansion. The open and filled circles 
represent fermion fields \ an d X, respectively, (a) In the LO, the time-like links (red arrows) in the 
fcrmion action Sp cancel with each other to leave a spatially local interaction, (b) In the NLO, the 
time-link links in Sf are canceled by the time-like links in a plaquette Sq (blue arrows) to leave a 
spatially non-local interaction. 



n \t(h 



c (3) I n=0 



x exp 



E 



/iV/2 \/ N/2 \ 

- 4 ( E Xs (x) X s {x + 4) J f J2 X.' {x + Z)Xs> (x) J 



\ E (V?(x)-Vf(x))+m.M{x) 



J=l,2 



(26) 



(27) 



where the factor in the second line fEa.([2l)])') consists of fermion interaction up to 2N- 
fermi terms. Hereafter we fix the flavor number N to 2 (monolayer case), so that the 
interaction is limited to 4-fcrmi term: 



4 0) = E 



\ E (V+(x)-V j -(x))+m.M(x) 
3=1,2 

-Y J M{x)M{x + l). 



(28) 
(29) 



C (3) 



Here the 4-fermi term in the second line fEa. (j29]) ) comes from the product of the temporal 
hopping terms V^{x) and V±(x) in Sf as shown in FigQJa). 

In the NLO (C^/? 1 )), two gauge formulations give some difference. With the compact 
formulation, the fermionic term Sf and one plaquette from Sg contribute to the NLO 
effective action 



4 1)c = ~ E E + 4) + (v; o v-)] 



(30) 



which come from the product of the temporal hopping terms V^~{x) and V 4 (x+j), and 
the plaquette U\{x)U±(x + j), as shown in FigJIJb) (and its Hermite conjugate). If we 
take the non-compact formulation, the NLO effective action reads 



s ^c = _ 2p{ _ 1) nJ2v i (x)-^J2 E M*)V4{*+3) + (const.), 

r(3) 3=1,2 
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(31) 



where the fermionic bilinear V4 is defined as 

V 4 (x) = x(x)x(x + 4) - x(x + 4) X (.t). (32) 
Here we have used the integration formulae 

— / dee = o (33) 

2"-7T J_ n7r 

/ d0 0e ±ie = T*(-l) n (34) 

2n?r J_ nv 

— / dO 6 2 e ±lB = 2(-l) n , (35) 
2n7r ,/_„„ 

where n is an arbitrary positive integer which specifies the interval of integration. We will 
show afterwards that the ultimate result with mean field approximation is independent 
of the choice of n. It should be noted that the effect of the gauge field in the £3 (or 
z)-direction starts to appear from 0(/3 3 ) in the strong coupling expansion. (In 0(/3 2 ), it 
only affects the constant term.) 

3.2. LO effective action with auxiliary field 

The non-local four-fermi interaction in the LO effective action (Eq.(T2§])) can be lin- 
earized by the extended Stratonovich-Hubbard (ESH) transformation [201 ] : 



e aAB 



= (const.) x J dipdip* exp[-a(|<y9| 2 - A<p - Bcp*)], (36) 



where A and B are fermion bilinears and a is a positive constant. Here we introduce a 
complex auxiliary field <fr(x) corresponding to the LO term, which transforms under the 
U(1) A rotation as <f>(x) — > e _2 ^A e ( x )0(x); thus we separate the real and imaginary parts 
of 4>(x) as 

4>(x) = <t>a(x) + ie(x)(f>T T (x). (37) 

The real part CT corresponds to the scalar operator M — xx> while the imaginary part 
cj>n to the pseudoscalar operator P — x^ e X- With the auxiliary field cj>, we can rewrite 
Eg. (1291) in terms of fermionic bilinears: 

5 eff= jEl^)! 2 - E x(x)V[<t>](x,y)x(y), (38) 

2,(3) s(3),y( 3 ) 

where the matrix element V[<f>](x,y) is defined as 

V[4>] (x, y) = y a>y [<f>(x) +4>*{x-i)- 4m,] - 1 £ Vj (x) \5 yx+] - S xy+] \ . (39) 

i=i,2 



This matrix element is written in the momentum space as 

(40) 



V[tft]{k,k') = \ [4>{p) + 4>*{~PV P - 4m»<J p ,o] + [(isinfci)<5 p)7r4 + (ismfc 2 )^ (4+1) 
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0.0 0.1 0.2 0.3 0.4 0.5 

1*1 

Figure 2: The free energy F e g((j>) in the lattice unit as a function of \<f>\ in the strong coupling limit 
(/3 = 0) and in the chiral limit (m = 0). 



where we denote p — k — k' . By integrating out the fermionic fields, the effective action 
is written only in terms of <f>: 



^° ) [0] = ^|^)| 2 -lndety[ 1 



(41) 



In order to diagonalize the matrix V[(j>] and calculate the determinant, here we per- 
form the mean-field approximation over <j) a and (/>„: 



4>ip) -> <PaSp,0 + ifaSp^fi 



(4+1+2) 



(42) 



in the momentum space. Thus the mass term in the effective action obtains the coefficient 
("effective mass") Mp = m* — </v/2, and V[<j>] yields the determinant 



lndetV^] = ^y^^ln [(T 1 ^)] , 



(43) 



with the effective bosonic propagator 
G(k; <f>) = 

and the two-dimensional momentum integration 



-i -l 



J'=l,2 



w dh r dk 2 

-rr 27T /__ 2TT 



(44) 



(45) 



From the above arguments, we can derive the free energy per unit cell (effective 
potential) analytically: 



ii k 



(46) 



Fig |2] illustrates the LO free energy F cff in the chiral limit (to = 0). From Ea. (l46l) . 
we find that F e g((f> —> oo) ~ |0| 2 due to the tree-level term, while F c s{4> — ^ 0) ~ 
const. + \4> | 2 In \4>\ 2 due to the fermion one-loop term. Therefore F^g takes the well-known 
"Mexican hat" shape, and dynamical chiral symmetry breaking always takes place in the 
strong coupling limit with the symmetry breaking pattern, U(l) v xU(l) A — > U(l) v . The 
potential minimum gives the chiral condensate 

\(xx)\=tr= | (0a) | = 0.240 (47) 

in lattice unit. 



3.3. NLO effective action 

There are two ways to linearize the NLO effective action Eqs. (l50|) and (f5T|) . One 
of them is perturbative expansion, which is rather straightforward; we use the fermion 
propagator given by the inverse matrix of V[<j)\. 



(0) _ 



(48) 



ik-(x-y) 



G(k; 



E Vj(y)(i sin kj 



where (• • -)(°) denotes the average over the LO effective action. Thus we obtain 

(V±(x))W=± f G(k^)sm 2 k 3 , (49) 



/ (1)C\ < ' ' / (l)NC\' ' 

which enables us to calculate (S^ ) and (S^ ) from Wick's theorem. The 
NLO term of the free energy in the compact formulation reads 



F 



(1)C 



ctr 



(0) 



E 

J=l,2 



G(k; (f>) sin 2 kj 



(50) 



On the other hand, we find by this process that the only term in S^ NG remaining over 
the average is 

"■E E [(Vf(x)Vr(x + 4) + (V+ o V-)] , (51) 

xO) 3=1,2 

and that all the other terms vanish up to 0{f3 1 ). Therefore the NLO effect in the non- 
compact formulation is twice larger than that in the compact formulation: 



t 
4 



F 



(1)NC 



(0) = 2F c W C (0) 



off 



(52) 



Another way to obtain the same results is to introduce an auxiliary field, which 
is rather indirect way than perturbative expansion, but more convenient for further 
applications, such as observation of collective excitations. By the ESH transformation, 
we introduce another complex auxiliary field X(x) corresponding to the NLO 4-fermi 
term V + V~ . Here we decompose X(x) into real and imaginary parts as A = Ai + i\2 and 
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perform mean-field approximation over them. The effective action in Ea. (|38[) is modified 
by the (compact) NLO term as 



w 2 + Ew 



j=l,2 



x(x)V[<P,\}(x,y) X (y), (53) 



where the modified matrix element V[4>, A] reads 



V[cf>,\](x,y) = -8 Xty [(j>(x) + 4>*(x-%)-lm* 



(54) 



-5(1 + ^) ^ Al/ \^-e-^ /2 ^+, 



i=i,2 



Thus the effective boson propagator in Eq. (|44|) is replaced by 



G(k;0,A) 



1 + — ) sm [ki- — 

3=1,2 



/3Ai 
2 



(55) 



and the effective potential in Ea. (j4"6l is modified correspondingly. 

By requiring the stationary condition dF^g + G / dX — 0, we can eliminate A by the 
relations 



G(k, <j>) sin ki cos h + 0(/3) = + 0(/3) 



A a = - / G(k 



I sm 



(56) 
(57) 



leading to the effective potential in the same form as the sum of Eas. (l46|) and (|50|) . As for 
the non-compact gauge formulation, the effective potential as the sum of Eqs. (l4"6")l and 
(fS"2"|) can be derived by the same process as for the compact formulation. Hereafter we 
mainly take the compact gauge formulation; non-compact results can be easily derived 
by substituting /3 in compact results by 2(3 up to the NLO terms. 

( 1 ) C NC 

Since the NLO correction term FlJ ' (4>) grows as \4>\ increases, the chiral conden- 
sate a is a decreasing function of /3. In other words, the chiral symmetry gets eventually 
restored as the coupling strength becomes weaker. Up to the linear terms in /3 and m, 
we can calculate a with the compact formulation as 

cr c (/3,m) ~ (0.240- 0.297^ + 0.0239 ma)a" 2 , (58) 

where we recover the lattice spacing a. If we take a -1 ~ a^l = 1.39 keV as a typical cut- 
off energy scale of this model, we obtain ct c (/3, m) ~ [(0.680 - 0.421/3 + ^r^) keV] 2 . 
Since the NLO term with the non-compact formulation is twice that with the compact 
formulation, the chiral condensate with the non-compact formulation er NC (/3) drops twice 
faster than a c {(i). The behavior of er c (/3) and <r NC (/3) around the strong coupling limit 
(/? = 0) is schematically shown in Fig0] The total fermion dynamical mass can be 
estimated from the coefficient of the effective mass term, 



v F aa 
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TO. 



(59) 
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Figure 3: The free energy F^a^ (<f>) up to NLO terms with the compact gauge formulation in the 
chiral limit. 




"0.072 



P 



Figure 4: Schematic picture of the behavior of exciton condensate <r(/3) in the strong coupling region. 
The compact (C) gauge formulation and the non-compact (NC) one coincides at /3 = 0, but tr(/3) with 
NC drops faster than that with C. Behavior in the region far from f) = (dashed curves) and the critical 
coupling value /3q cannot be examined by strong coupling expansion, so that here we employ the results 
suggested by lattice Monte Carlo method [l7j| . 



If we take a ~ a Hc and the compact formulation, we can estimate Mp ~ (0.523 — 
0.623/3) eV + 3.05m. 

Since strong coupling expansion is valid only for small-/? region, it is no longer possible 
to investigate the nature of the chiral (scmimctal-insulator) transition: the critical cou- 
pling value /3c , the order of the transition, and the critical exponents for the transition, 
unless some resummation techniques are introduced. Monte Carlo simulations provide 
some clues to such critical behavior [l7j : with the same lattice model employed in our 
studies, the phase transition is estimated to be of second order with flc = 0.072 ± 0.003 
with the non-compact formulation. With the compact formulation, it is estimated to be 
of first order with (i c = 0.42 ± 0.01. 
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Due to the finite-size effect of lattice, lattice Monte Carlo simulations cannot reach 
the strong coupling and chiral limit. Here we approximately compare our results with 
the Monte Carlo results in the strong coupling and chiral limit, by extrapolating the 
Monte Carlo results with the equation of state which has been conjectured in analogy 
with the (3+l)-dimensional QED [I?}. The chiral condensate reaches the value around 
0.24 (in lattice unit) in the strong coupling limit by Monte Carlo simulation, both with 
the compact formulation and with the non-compact one. By increasing the value of /?, it 
can be seen that cr (/3) drops faster than cr c (/3). These results are consistent with the 
results obtained from our work. We expect that our analytical studies, applicable around 
the strong coupling and chiral limit, and the Monte Carlo simulation results, applicable 
within the finite coupling strength and finite mass region, are complementary to each 
other. 



4. Collective Excitations 



Here we consider two kinds of collective excitations, or the fluctuations of the order 
parameter <p(x) around the symmetry broken state (</>) = — cr: the phase fluctuation mode 
corresponding to ^(x) and the amplitude fluctuation mode corresponding to 4> a (%)-> 
which we name "7r-exciton" and "cr-exciton" respectively in analogy to the pion and the 
cr-meson in QCD. Propagator of the a-exciton mode (a = cr, n) can be derived from the 
effective action 5 e ff[<^] (with the auxiliary field method) as 



Is 



(50 Q (x)50 Q (y) 
+ Tr 



r =-cr,<4„=0 



t {x) d<j) a {y) 



where we have used the formula 
d 2 \ndetV 

— Ir 







+ Tr 


d(f>id4>2 



(60) 
(61) 

(62) 



and the trace runs in the (2+l)-dimensional space-time. The second term in Eq. (|6"2"|) 
vanishes because the matrix element V[<^](x, y) is linear in </>. The second term in Ea. (|61[) 
corresponds to the fermion one-loop diagram with two </> a -legs, as shown in FigfSJ Thus 
the propagators can be written in the momentum space as 



_ i 1 1 + cosh c 
Dj (p, = 



H(k, P] a)G(k-a)G(k + p;a), 



(63) 



where G is the "effective" propagator defined in Ea. (j4~4"|) (in LO) and Eg. (1551) (in NLO), 
with <f> and A substituted by their expectation values determined by stationary conditions. 
The effect of two vertices in the diagram is represented by 



i?(k, p;cr) = ± (m* + '- 



j=l,2 
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where the ± sign corresponds to a = tt and a respectively. The dispersion relation for it 
and <7 is derived from the pole of the propagator: by restoring the temporal scaling by 
v F , we obtain the dispersion relation 



(65) 



Specifically the mass, or the excitation energy, of these modes are given by the energy 
at zero momentum: M a ^ = w ff!7r (p = 0). 




Figure 5: The excitonic propagator consists of a fermion one-loop diagram with two </v,7r legs, 

with the overall momentum p. 



n-exciton mode 

By considering the pole equation for 7r-exciton D^ 1 (p = 0,u;J), the mass M w 
v F u)*(p = 0) is given by 



coshw! 1 



G(k;a) 



1. 



(66) 



With the help of the gap equation [dF e s((j)) / = 0, the integration part is given 
as 

M F (m = 0) 



G(k;a)=2- 
k M F (m) 



(67) 



where Mp(m = 0) is dynamically generated mass in the chiral limit Mp(m = 0) = v F a/2. 
Thus Ea. drJBl yields coshojj = 1 + 2m/M F (m = 0). By taking the first order in the 
fermion bare mass m and restoring the lattice unit a, the mass M 7r (m) of this mode 
reads 



Mir 



a V M F (m = 0) 



(68) 



Thus the 7r-exciton serves as the massless Nambu-Goldstone (NG) boson related to 
the spontaneous breaking of the U(1)a symmetry, which is emergent in the low-energy 
approximation. As long as < m < 2 meV is satisfied, < M F holds, so that the 
7r-exciton serves as the lightest mode in the system. The 7r-exciton behavior derived 
here is reliable as long as the excitation energy M n is within the scale of the low-energy 
approximation, or, in other words, m <C M F (m = 0). 
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The relation M„ oc y/m is similar to the Gell-Mann-Oakes-Renner (GMOR) relation 
for the pion obtained from current algebra in QCD 22, Hj]. Considering the infinitesimal 
local U(1)a transformation 

S\(x) = ie(x)a(x)x(x) (69) 
and taking the transformation of (P(x)) in the present system, we obtain 

= 5(P(x)) = | j[d x d x ][d8]5[P{y)e~ s ] (70) 
= (SP(y) - P(y)5S) , 
where the infinitesimal transformations of the pseudoscalar density and the action read 
SP(y) = -2a(y)M(y) (71) 

ss = Y, a ( x )[- d n J f al ( x ) + 2m * p ( x )]> ( 72 ) 

x (3) 

with the axial current defined as J axial (a;) = ^e{x)(V^ (x) — V+{x)). By applying 
the functional derivative 5/5a{x) to both sides of Eg. ([70]) . we obtain the axial Ward- 
Takahashi (WT) identity 

(P(y) (d„jf*\x) - 2m*P(x)) - 2M{y)5 xy ) = 0. (73) 

Saturating this WT identity by the pole of the 7r-exciton, we obtain 

p -ipx 

dM0\P(y)\Ap)) — (Ap)\P^\x)+2m*P(x)\0) - -2aS xy (74) 

If we take x = y = and m* = 0, the matrix element 

(0\P(0)fr( p )) = --jL + O(m) (75) 

is given in the leading order of m, where the temporal "pion decay constant" F£ is 
defined by the matrix element, (0| J| xlal (0)|7r(p)) = 2F£p^. On the other hand, if we 
take x = y and finite to*, we obtain 

2MlFl + 2m,(0|P(0)K(p)) = 0. (76) 

Thus, from Eas. (175p and (|76[) . we obtain, in the leading order of to, 

{FlM^f = ma, F :=^=a-^, (77) 

where a takes the value in the chiral limit (to = 0), and the lattice unit a is again 
restored here. The first equation in Eq. ([77|) , which tells that the mass of 7r-exciton is 
proportional to the square root of the fermion bare mass, is indeed the same form as the 
GMOR relation. 
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4-2. a-exciton mode 

As for the cr-exciton, we obtain its mass 



(1.30 - 0.47/3)- 



22.6m 



(78) 



by solving DT^(0,iM a /v F ) = numerically. This reduces to M a ~ (5.47 - 1.97/3) eV + 
22.6m when we take a ~ a Hc . This value is comparable to the cutoff energy scale of 
the present lattice, Ea = v F Tr/a — 13 eV. Although the cr-exciton does not have width 
in the present one-loop approximation, it would eventually decay into a 7r-exciton pair 
in higher orders. This is analogous to the situation for the broad er-meson in QCD 
However, since this value of M a is beyond the scale of low-energy approximation, its 
numerical accuracy is not reliable in the present approach. In order to be more precise, 
the model without the low-energy approximation, such as one preserving the structure 
of the original honeycomb lattice, is required. Shown in FigjS] is an illustration of the 
spectrums of the fermion and collective exciations obtained in this study. 



•E A /2 



E A /20 







Mo 



M F 



(3=0 
m=0 



0<3«1 
m=0 



0<p«1 
m~E A /10 3 



Figure 6: A schematic picuture of the fermion excitation energy Mp, the 7r-exciton mass M v and the 
cr-exciton mass M a obtained from the strong coupling expansion. = and m = correspond to the 
strong-coupling limit and chiral limit, respectively. The ultraviolet cutoff for the energy, Ea, is given by 
v F TT/a. 



5. Conclusion 

In this paper, we have investigated the behavior of the chiral symmetry of monolayer 
graphene in the strong coupling regime analytically by using the low-energy effective 
field ("reduced QED") model of monolayer graphene. A detailed explanation about the 
strong coupling expansion of the lattice-regularized model, which was first applied for 
graphene model in our previous paper (2lj . is given. "Chiral symmetry" (sublattice on- 
site energy balance) of the original honeycomb lattice, of the continuum effective theory, 

18 



and of the staggered fermion formulation, are compared. Since the compact U(l) gauge 
theory suffers from anomalous monopole condensation, we compared the results with 
compact and non-compact formulations of the gauge field action. As a result, whichever 
formulation we take, chiral symmetry is spontaneously broken in the strong coupling limit 
with equal value of chiral condensate. As the coupling strength becomes weaker, chiral 
condensate with the non-compact formulation drops faster than that with the compact 
formulation. These results up to NLO terms in strong coupling expansion agree with the 
Monte Carlo simulation results, extrapolated to the strong coupling limit by the equation 
of state. 

We have also estimated the energy scale of the emergent fermionic and collective 
(bosonic) excitations, and have found that the fermionic quasiparticle and the 7r-exciton 
are within the scale of the low-energy approximation, as long as the on-site energy dif- 
ference between two sublattices is within that scale. On the other hand, the mass of 
cr-exciton is quite large even in the chiral limit, and it is beyond the low-energy approx- 
imation. Thus one can say that the cr-exciton needs more investigation with a model 
without low-energy approximation, such as one preserving the structure of the original 
honeycomb lattice. Moreover, since the physical behavior of multilayer graphene, in 
which interlayer hopping amplitude depends on the stacking pattern between layers, is 
also of great importance, an analytic investigation with such a honeycomb lattice model 
is now strongly required. We expect that the U(l) gauge theory on the honeycomb lattice 
can be suitable both for analytic calculations and for numerical simulations. 

Extension of this analysis to the finite temperature and finite density region is of 
another important question. In order to extend our analysis to the finite temperature 
region, we have to take a temporally finite-size lattice and incorporate fluctuation effects 
beyond the mean-field approximation. This procedure would also enable us to treat the 
dynamics in the temporal direction distinctively so as to investigate the renormalization 
effect on the Fermi velocity. Finite temperature and chemical potential analysis would 
be required not only for the chiral transition but also for transport properties, such as 
electric conductivity and Hall effect. 
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